Transmission probability through a Levy glass and comparison with a Levy walk 
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Recent experiments on the propagation of light over a distance L through a random packing of 
spheres with a power law distribution of radii (a socalled Levy glass) have found that the transmission 
probability T cx 1/L^ scales superdiffusively (7 < 1). The data has been interpreted in terms of 
a Levy walk. We present computer simulations to demonstrate that diffusive scaling (7 ~ 1) can 
coexist with a divergent second moment of the step size distribution (p(s) oc 1/s^^" with a < 2). 
This finding is in accord with analytical predictions for the effect of step size correlations, but 
deviates from what one would expect for a Levy walk of independent steps. 

PACS numbers: 05.40.Fb, 05.60.Cd, 42.25.Bs, 42.25.Dd 
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I. INTRODUCTION 

A random walk with a step size distribution that has 
a divergent second moment is called a Levy walk [ll-Q- 
A Levy glass is a random medium where the separation 
between two scattering events has a divergent second mo- 
ment. The term was coined by Barthelemy, Bertolotti, 
and Wiersma for a random packing of polydisperse 
glass spheres. They measured the fraction T of the light 
intensity transmitted through such a random medium 
in a slab of thickness L, and found a power law scal- 
ing T cx with a superdiffusive exponent 7 ~ 0.5 

— intermediate between the values for ballistic motion 
(7 = 0) and regular diffusion (7 = 1). 

The simplest theoretical description of propagation 
through a Levy glass neglects correlations between sub- 
sequent scattering events. The ray optics of the problem 
is then described by a Levy walk, with a power law step 
size distribution p{s) cx l/s^'^°', < a < 2. The experi- 
ment [3] was interpreted in these terms, with a = 1 and 
7 = a/2 the expected transmission exponent. 

Correlations between scattering events in a Levy glass 
dominate the dynamics in one dimension 0,0. Ahhough 
correlations were expected to become less significant with 
increasing dimensionality 0,111, Buonsante, Burioni, and 
Vezzani 9] have calculated that the transmission expo- 
nent 7 should remain much larger than would follow from 
a Levy walk with uncorrelated steps. In particular, a sat- 
uration at the diffusive value 7 = 1 for a > 1 is predicted 

— even though the second moment of the step size dis- 
tribution becomes finite only for a > 2. 

To test these analytical predictions for the effect of 
correlations, we have simulated the transmission of clas- 
sical particles through a Levy glass, confined to a slab 
of thickness L. Both a two-dimensional (2D) system of 
discs is considered and a three-dimensional (3D) system 
of spheres. Wc find a power law scaling T{L) cx 1/L'' 
with an exponent 7 that lies well above the ^ — a/2 line 
expected for a Levy walk. In particular, we obtain a sat- 
uration of 7 at the diffusive value of unity well before the 
a = 2 threshold is reached of a divergent second moment. 



The outline of the paper is as follows. Since our aim 
is to compare the Levy glass simulations with the pre- 
dictions for a Levy walk, we need analytical results for 
uncorrelated step sizes. These are summarized in the 
Appendix and referred to in the main text. We start 
off in Sec. |lT] with a description of the way in which we 
construct and simulate a Levy glass on a computer. The 
results presented in that section are for 2D, where the 
largest systems can be studied. We turn to the 3D case 
in Sec. IIIII and compare with the experiments We 
conclude in Sec. HVl 



II. LEVY GLASS VERSUS LEVY WALK 

A. Construction 

A Levy glass [3,IToi is a random packing of transparent 
spheres with a power law distribution of radii. 



n{r) cx l/r 



(2.1) 



Light propagates without scattering (ballistically) 
through the spheres and diffusively (mean free path ^mfp) 
in the region between the spheres. The probability to en- 
ter a d-dimensional sphere of radius between r and r + dr 
is proportional to the fraction n{r)dr of spheres in that 
size range, multiplied by the area cx r*^"^. The ballistic 
segments (steps) of a ray inside a sphere of radius r have 
length s of order r. The sphere radius distribution (12.11) 
therefore corresponds to the step size distribution [ll| 



p{s) cx l/s""" 



with /? = a + d - 1. 



(2.2) 



Particles propagating through a Levy glass therefore 
have the same distribution of single step sizes as in a 
Levy walk, but the joint distribution of multiple step 
sizes is different: While in a Levy walk the steps are all 
uncorrelated (annealed disorder), in the Levy glass the 
configuration of spheres is fixed so subsequent steps are 
correlated (quenched disorder). 
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FIG. 1: Two-dimensional Levy glass, consisting of a random 
packing of discs with a power law distribution of radii (a = 
0.7, / — 0.86, rmax/'"min ~ 100). The blue region defines 
a slab of thickness L. This is the unconstrained geometry, 
because the maximum disc size can be larger than L. 



We discuss in some details the construction of the 2D 
Levy glass, see Fig. [1] — the 3D version is entirely anal- 
ogous. We start by generating discs of (dimensionless) 
radius 



A; = 0, l,2,...fc: 



k 



1 + — (r^ - 1) 

h. V' max 



(2.3) 



The /cmax + 1 discs have radii ranging from Tmin = 1 to 
Tmux S> 1, and in this size range their distribution follows 
the power law (|2.ip . The average area of a disc is 



(A) 



|2-/?| 



max(l, 



(2.4) 



The entire Levy glass occupies an area of dimension 
X in the x — y plane, with periodic boundary con- 
ditions and W about 10-100 times larger than Tmax- For 
a random packing we place the discs at randomly chosen 
positions in the order k = 0,1,2,... (so starting from 
the largest disc). If disc number k overlaps with any 
of the discs already in place, another random position 
is attempted. For each disc some 10'' attempted place- 
ments are made. If they are all unsuccessful, the entire 
construction is started over with a smaller value of fcmax- 
The density of the packing is quantified by the filling 
fraction 



/ = fc„ax(A)/W^2 



(2.5) 



For each simulation we strove for maximal /, by maxi- 
mizing fcmax- The maximal filling fraction increases with 
increasing ratio rmax/''min, as illustrated in Fig. [5J For 
the smallest a, below about 0.4, we could not reach as 
dense a packing as for larger a, basically because there 
are too few small discs. Somewhat larger filling fractions 
would be reachable by moving the discs after placement, 
but we did not attempt that. 
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FIG. 2: Filling fraction of the 2D Levy glass as a function 
of the ratio rmax/j'min of largest and smallest disc size, for 
several values of the parameter a. 




FIG. 3: Step size distribution for a random packing of discs 
with radius distribution pTTj) (for 13 = 2.2, so a = 1.2). The 
numerical results are shown for two values of the maximum 
disc radius (rmax/rmin ~ 10^ and 10^, with / = 0.83 and 0.80, 
respectively) . The black solid line is the expected distribution 



B. Dynamics 

The ballistic dynamics inside the spheres consists of 
chords of varying length s traversed in a time s/v. The 
diffusive dynamics in between the spheres is modeled by 
a Poisson process: isotropic scattering in a time inter- 
val dt with probability wdt/Zmfp. The mean free path 
^mfp = ^"111111/2 is chosen such that there is, on average, one 
scattering event between leaving and entering a sphere. 
We take the same refractive index (and velocity v) inside 
and outside the spheres, so the ray is not refracted at the 
interface. 

In Fig. [3] we show the step size distribution p{s) for 
a 2D Levy glass with disc radius distribution (|2.ip . for 
/3 = 2.2. It follows closely the Levy distribution ()2.2p . 
with the expected parameter value a = /3— 1 = 1.2 (solid 
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FIG. 5: Diffusion coefficient (12.61) in the Brownian regime, 
estimated from tlie large-t slope of the mean square displace- 
ment (corresponding to the large-t saturation value in Fig. 
|4]). Each set of colored data points represents one value of 
Levy glass with different values of rmax/»'min, at fixed a — 1.2. a, with different values of rmax/rmin. The power law scaling 

(|2.7p (red dotted lines) determines the scaling exponent 70. 



FIG. 4; Time dependence of the mean square displacement 
(divided by t so that saturation indicates diffusive scaling). 
The curves are the results of a numerical simulation in a 2D 



line) . 

We do not find the pronounced oscillations in p(s) 
which in Ref. [13 complicated the determination of a. 
These oscillations appear due to coarse graining of the 
disc size distribution n(r) and vanish if a finer distribu- 
tion of disc sizes is used. 

The time dependence of the mean squared displace- 
ment (Ar(t)^) is shown in Fig. |4l for the same a = 1.2. 
A particle was started at a random position r(0) in the 
inter-disc region, and then its position r(i) at time t 
(either inside or outside a disc) gives the displacement 
Ar(i) = |r(t) — t(0)|. The average (•••) is over some 
10'* initial positions. In accord with previous simulations 
[1, [1^, regular (Brownian) diffusion with (Ar(t)^) cx t is 
reached for times t ^ r^^-y^jv = ^d, set by the time needed 
to traverse the largest disc. For t < tu the mean squared 
displacement increases more rapidly than linearly (su- 
perdiffusion) . 

The limiting slope of the mean square displacement 
for t 3> gives the diffusion constant in the Brownian 
regime. 



t-s-oo 2at 



(2.6) 



As shown in Fig. [51 this diffusion constant has a power 
law dependence on Tmax, 



1-7D 
max 1 



(2.7) 



with < 7£) < 1. (For the smallest a — 0.2 no clear 
power law scaling was observed.) 



geometry. One way is to obtain the slab from the entire 
Levy glass by cutting out the region < x < L (blue 
strip in Fig. [T]). We call this an unconstrained geome- 



try, because 



is not constrained to be smaller than 



L. The alternative constrained geometry (used in the ex- 
periments [3]) requires that the spheres all lie fully inside 
the slab, thereby restricting Tmax < L/2. We consider 
the transmission probabilities in the unconstrained and 
constrained geometries in separate subsections, both for 
2D. (Results for 3D are presented in the next section.) 



D. Unconstrained geometry 



A 



lower limit Tbaii to the transmission probability 
Tuncon in the unconstrained geometry follows by consid- 
ering only ballistic rays, which pass through the region 
< a; < L without a single scattering event. As ex- 
plained in the Appendix, see Eq. (|A7|) . this probability 
is directly related to the step size distribution. 



Jball — 



1 



dx 



dsp{s). 



(2.8) 



We take the step size distribution p.2|) with an upper 
cutoff at Smax — ''max ^ L and a lower cutoff at Smin — 1. 
Then Eq. (IS^Hl) evaluates to 



ball 



C. Transmission probability 

For the transmission problem we need a slab of variable 
thickness L. We distinguish two ways of constructing this 



1 



for < a < 1, 
for 1 < a < 2. 



(2.9) 



Since Tum < T. 



< 1 we can immediately conclude 



that r„ 



1 for < a < 1. For 1 < a < 2 the 
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FIG. 6: Transmission probability Tuncon through a 2D un- 
constrained Levy glass, for different values of the maximum 
disc radius rmax. The dotted line is the predicted scaling 
Tuncon oc L^'"' in the r^ax -> oo limit. 



power law scaling Tuncon oc 1/ L'^ must satisfy j < a — 1. 
This holds irrespective of correlations between multiple 
steps, since these cannot affect Tbaii- If we neglect these 
correlations, we may equate Tu„con to the transmission 
probability T^q of a Levy walk with equilibrium initial 
conditions (see App. I A 3")) . In view of Eq. (|A9I) . this leads 
to 7 = a — 1. We believe this result to be quite robust, 
since even if correlations do play a role, it is likely that 
they slow down the superdiffusion ,7, 8], so they would 
not lead to a smaller 7. 

In Fig. [6] we show the L-dependence of Tuncon for two 
values of a, resulting from a numerical simulation of an 
unconstrained 2D Levy glass. This is data up to rmax = 
10'' for a = 1.1 and up to r^ax = 10"^ for a = 1.5, which 
is at the upper limit of our computational resources. As 
expected from the Levy walk (Fig. [T5|) . the convergence 
to the Tmax 00 limit is very slow, and we are not able 
to conclusively test the predicted asymptote. 



E. Constrained geometry 

For the construction of a constrained Levy glass we 
limited the maximum disc radius to r„iax = T/4 and 
ensured that all discs fit inside the slab of thickness L. 
The corresponding random walk would be a truncated 
Levy walk with maximum step size s„iax — L/2. From 
the analysis in App. IA4bl we would therefore expect a 
T cx 1/L"/^ scaling of the transmission probability — if 
correlations between step sizes would not matter. 

In Fig. [7] we show the scaling of the transmission prob- 
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FIG. 7: Transmission probability through a 2D constrained 
Levy glass as a function of the thickness of the slab, for dif- 
ferent values of the step size exponent a. The dotted lines are 
a linear fit to the data points, determining the transmission 
scaling exponent 7. ( The data is split over two panels, to 
avoid overlap.) 
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FIG. 8: Exponents 7 and 70, governing the scaling of the 
transmission probability (|2.10p (crosses) and diffusion con- 
stant (|2.7|l (circles) . These are the results of a simulation of a 
2D constrained Levy glass (see Figs.[5]and[7]). The red dashed 
line is the prediction (|A20[) for a Levy walk with nonequilib- 
rium initial conditions. 



ability. 



Tex l/L'' 



(2.10) 



as it follows from the simulation. The power law scaling 
applies to somewhat less than two decades in L for a ^ 
0.5 (lower panel), and to one decade for smaller a (upper 
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FIG. 9: Comparison of the a-dependence of the transmission 
exponent 7 for a 2D and 3D Levy glass. Both data sets lie 
well above the 7 = 0/2 line of a Levy walk. 

panel). In Fig. |8]we give the resulting exponent 7 as a 
function of a. 

In the same figure we show the scaling of the diffusion 
exponent 7D, from Eq. (12. 7|) . (There we could only ob- 
tain a power law scaling for a ^ 0.4.) As expected from 
the identification of T ~ D{L)/L cx one has in 

good approximation 



7 = Id- 



(2.11) 



III. COMPARISON WITH EXPERIMENTS 

The numerical data shown so far was for a 2D Levy 
glass of discs. We have also performed simulations for 
a 3D Levy glass of spheres, in the constrained geometry 
with Tinax = L / A. We went up to L/r,„in — 1132 for a < 
0.8 and up to L/r^i^ — 800 for a = 1 and 1.2. (Larger 
values of a could not be simulated reliably.) Although 
the systems are smaller in 3D than in 2D, the results 
are quite similar, see the comparison in Fig. [3] of the a- 
dependence of the transmission exponent 7 for a 2D and 
a 3D Levy glass. In particular, for both 2D and 3D the 
results for 7 lie well above the ^ — a/2 line. 

We can now compare directly with the 3D experiments 
which obtained 7 = 0.5 within experimental accuracy 
for a — 1. Our simulation, in contrast, gives for a ^ 1 
a value for 7 which is about 50% higher. We cannot 
attribute the difference to finite-size effects, since the 3D 
simulation reaches the same range of system sizes as the 
experiment. There are aspects of the experiment which 
are not present in the simulation (notably absorption), 
but we believe that the difference is mainly due to an 
irregularity in the experimental sphere size distribution. 

To visualize the irregularity we plot in Fig. [TU] the 
quantity 
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FIG. 10: Sphere volume distribution used in the experiment 
U (red solid histogram) and for an a = 1 Levy distribution 
(green dashed histogram). 



which is the cumulative volume enclosed by spheres with 
radii greater than r. This is a decreasing function of r, 
from l^(rmin) = Vo (the total sphere volume) down to 
V^(rniax) — 0. For the Levy distribution with a = 1 in 
3D we have n{r) oc r~^, cf. Eqs. (|2.ip and (|2.2I) . hence 
V{r) should decrease linearly as a function of logr. 



V{r) 



Vo 



l0g(r„iax/rmin) 



log(r/rmax). 



(3.2) 



V{r) = Itt / r'^n{r')dr', 



(3.1) 



As shown in Fig. llOi the experimental sphere size dis- 
tribution differs markedly from the expected Levy form 
p.2p . Rather than a single linear dependence of V{r) on 
log r, there are two piecewise linear dependencies with a 
different slope, joined with a kink at r « 50 /im. This 
irregularity has the effect of reducing the transmission 
exponent 7, essentially by mimicking a system with a 
smaller value of a. 

To demonstrate the effect of the kink on the trans- 
mission exponent we have simulated the experiment by 
constructing a random packing of spheres with the exper- 
imental size distribution (red solid histogram in Fig. [T0| . 
All spheres were constrained to fit inside a slab of thick- 
ness L. (We took r,„ax = L/2.1 for this simulation.) We 
found 7 — 0.57. If instead we used the proper Levy size 
distribution (green dashed histogram), keeping all other 
parameters the same, we found 7 — 0.72. We believe this 
resolves the issue. 



IV. CONCLUSION 

In conclusion, we have found that the superdiffusive 
scaling T cx 1/ L'^ of the transmission probability through 
a Levy glass, constrained to a slab of thickness L, devi- 
ates substantially from what one would expect for a Levy 
walk. Most significantly, the diffusive scaling (7^1) can 
coexist with a divergent second moment of the step size 
distribution (a < 2). 



6 



As a consistency check on our simulations, we have also 
calculated the diffusion constant D from the long-time 
limit of the mean-square-displacement in an unbounded 
Levy glass, as a function of the maximum disc size rmax- 
We find £'(7'max) c>c r^^g^" , with jjy ~ 7, as expected for 
a diffusive transmission probability T ~ D{L)/L with a 
scale dependent diffusion constant. 

Qualitatively, our finding that diffusive scaling of T 
can coexist with a divergent second moment of p(s) is 
consistent with analytical calculations for d = 1 [51 and 
d = 2, 3 0. Quantitatively, we are not in agreement: Ref. 

finds that 7 increases monotonically for d = 2 from 
7 = at a = to 7 = 1 for a > 1, while our simulation 
gives a nonmonotonic a-dependence of 7, with a satu- 
ration for a > 1.5 (see Fig. [5]). The system considered 
in Ref. f9| is quasiperiodic (a Levy quasicrystal) , rather 
than the random Levy glass studied here. Further study 
is needed to see whether this difference is at the origin 
of the different transmission scaling, or whether the dif- 
ference is due to a very slow convergence to the infinite 
system-size limit (which we consider more likely). 
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Appendix A: Transmission probability of a Levy 
walk 



1. Formulation of the problem 

We consider a random walk along the a;- axis with the 
power law step size distribution 

p(s) = -(-)'^%(s-so). (Al) 

Sq \ S / 

(The function 9{s— sq) equals 1 if s > sq and if s < sq.) 
Subsequent steps are -f s or — s with equal probability and 
independently distributed. The probability density p{s) 
decays as l/s^+" with a > 0, starting from a minimal 
step size sq > 0. In between two scattering events the 
walker has a constant velocity of magnitude v. This ran- 
dom walk is called Brownian or diffusive for a > 2, Levy 
Q or superdiffusive for 1 < a < 2 and quasiballistic for 
< a < 1. 

The walker enters the segment < a; < L by passing 
through a; = at time U and then stays in that segment 
until time tf. If at it exits through x = L we say the 
walker has been transmitted through the segment. We 
seek the dependence of the transmission probability T on 
the length L of the segment, for L l^. For a Brownian 
walk, the scaling is inverse linear: T cx 1/L if a > 2. 
For a Levy walk we expect a slower power law decay, 
T cx 1 /L^ with 7 < 1 . The question is how 7 varies with 
a < 2. 




FIG. 11: Trajectories x(t) of a random walk, with scattering 
events indicated by red dots. All trajectories enter the seg- 
ment < a: < L (between dotted lines) at a; = 0. Trajectories 
a,b,a',b' are transmitted through x = L, while trajectories 
c, c' are reflected through a; = 0. The transmission probablity 
Tcq averages over all trajectories (equilibrium initial condi- 
tions), while Tiioneq averages only over trajectories such as 
a',b',c' that have a scattering event upon entering the seg- 
ment at a; = (nonequilibrium initial conditions). 



The answer depends on how the walker is started off 
initially. Following Barkai, Fleurov, and Klafter [T^l, we 
distinguish equilibrium from nonequilibrium initial con- 
ditions. fSee Fig. 1111 ) For equilibrium initial conditions, 
the walker starts off from x — —00, so that it crosses 
a; = at some random time between two scattering 
events. For nonequilibrium initial conditions, the walker 
starts off from a; = with a scattering event. We denote 
the transmission probabilities in these two cases by Teq 
and Tnoneq, respectively, and consider the two cases in 
separate subsections. 



2. Nonequilibrium initial conditions 



The transmission probability Tnonoq from a; = to a; = 
L for a Levy walk that starts off with a scattering event 
at X = has been calculated by several authors |13l- 15|. 
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We give the most general solution of Buldyrev et al. 

They assume that the walker starts with a scattering 
event at an arbitrary point Xi in the segment (0, L) and 
calculate the probability P{xi) that the walker exits the 
segment through x = L. For L ^ sq and Xi ^ sq their 
solution [l5| can be written in the compact form 



B(.x,/L,a/2,a/2) 



B{l,a/2,a/2) 
in terms of the incomplete beta function 



B{x, a, b) 



y 



(A2) 



(A3) 



Since B{x, a, b) — x'^/a for a; — ^ 0, one arrives at the scal- 
ing Tnoneq oc first obtained by Davis and Marshak 
from basic considerations uM. 
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The prefactor of the power law scahng cannot be ob- 
tained directly from the solution (|A2p . because of the 
limitation that Xi 3> sq. For < a < 1 we can work 
around this limitation by considering the first step sep- 
arately. The walker starts off at a; = with a step to 
xi > 0, chosen randomly from the distribution ()A1|) of a 
Levy walk. If xi > L the walker is transmitted with unit 
probability. Otherwise, it is transmitted with probability 

We thus can calculate Tnoneq from 



dxip{xi)+ / dxip{xi)P{xi). (A4) 



For a < 1 the mean step size diverges, so the region 
Xi < So is insignificant and we can use Eq. (|A2[) for 
P{xi). The result is 



noncq 



B{so/L,a/2,l + a/2) 
a/2,1-1- a/2) 
L»so /so\"/2 4r(a) 

\l) ar2(a/2)' 



(A5) 



While the exponent a/ 2 holds for any < a < 2, the 
prefactor is accurate only for < a < 1. (For a > 1 we 
would need to know P{xi) within the region xi ^ sq in 
order to calculate the prefactor.) 



because the walker is more likely to cross x = during 
a long step than during a short step, so long steps carry 
more weight in q{x) than they do in p(s). Indeed, for 
1 < a < 2 the first moment of q{x) is infinite while the 
first moment oi p{s) is finite. 

Substitution of Eqs. (|X2|) and into Eq. gives, 
for L > So, 



7rr(a) 



asin(a7r/2)r2(a/2)' 



for 1 < a < 2. 



(A9) 

This scaling T^q oc 1/L" ^ holds in the superdiffusive 
regime 1 < a < 2. In the quasiballistic regime the first 
scattering event is at xi > L with unit probability. 



T,, 



1, for < a < 1. 



(AlO) 



The value a = 2 at the border between a Brownian 
walk and a Levy walk requires separate consideration. 
While Tnoneq OC 1/L for a = 2, the transmission proba- 
bility (|A6p has a logarithmic enhancement. 



L 



1 



2 so 



for a 



(All) 



A similar but different scaling cx L ^Vln L has been as- 
sociated with the a = 2 Levy walk in Ref. [l^ 



3. Equilibrium initial conditions 

For equilibrium initial conditions the walker crosses 
a; = at a random time between scattering events. The 
first subsequent scattering event is at a point xi > 0, 
with probability density q{xi). li xi > L the walker is 
transmitted with unit probability, if < xi < L the 
transmission probability is P{xi). Hence 



dxiq{xi)+ / dxi q{xi)P{xi). (A6) 



The probability density q{x) is determined from the 
step size distribution. 



1 f°° 
(s) Jx 



(A7) 



This relation between the distribution p{s) of the dis- 
tance s between subsequent scattering events and the dis- 
tribution q(x) of the distance x from an arbitrary point 
to the next scattering event holds for any random walk 
with a finite average step size (s) = sp{s) ds. For the 
step size distribution (jAip one has 



q{x) 



1 



So 



aso \max(x. So) 



for a > 1. 



(A8) 



As emphasised in Ref.[i3, the distribution q{x) oc l/x" 
decays more slowly than the distribution p{s) oc l/s^"*"" 



4. Truncated Levy walk 



A truncated Levy walk has step size distribution 

l+a 

So \ s ' 



a /So\i+" 

Ptrunc(s) = — — 9{s- So)6'(s„iax - s), (A12) 

Sn \ s / 



with a maximum step size Smax 3> so- The root- mean- 
squared displacement a after a single step then has a 
finite value. 



a l-a/2 a/2 
'-'max "^0 ' 



2-a 



(A13) 



much smaller than Smax for a < 2. 

The transition from a truncated Levy walk to a Brow- 
nian walk requires ngtops ^ 1 of steps, given by [l^, [l3| 

(2 - a)3 

r^stcps - ^ ^(Smax/So)"- (A14) 

a 

The corresponding root-mean-squared displacement 



cps 



(2 — a)smax is of order Smax for all a < 2. We 
conclude that we have regular (Brownian) diffusion over 
a distance L if Smax ^ L. 

The transmission probability P{x) for a walker start- 
ing with a scattering event at a point x inside a slab of 
thickness L (further than Smax from the boundaries) thus 
follows the usual diffusive scaling. 



P{x)=x/L, if x,L — x^ 



(A15) 
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a. Equilibrium initial conditions 

For equilibrium initial conditions the distribution q{x) 
of the first scattering event follows from Eq. (jA7[) . with 
p replaced by Ptrunc- Substitution into Eq. (jA6p then 
determines the transmission probability (for L > Smax), 



T — 

cq — 



dx q{x)P{x). 



(A16) 



Eq. (IA15[) gives P{x) only for x ^ Smax- We will use this 
expression also for x < Smax, and then test the approxi- 
mation by comparing with numerical simulations in Sec. 

If we substitute P{x) = x/L wc find 



T -^^ 

cq — - ^ 



a s„ 



max 



2i 2 — a s„ 



for < a < 1 or 1 < a < 2. For a = 1 or a 
are logarithmic factors, 



^^'^ 2Lln(s„ 



T,. 



So Sn 

2L 



ax /so) 
lll(Smax/So) 

max 



for a = 1, 



for a 



(A17) 
2 there 

(AlSa) 
(A18b) 



For fixed Smax the diffusive 1/L scaling holds. An 
anomalous scaling appears if the maximum step size 
Smax = cL is a fixed fraction c < 1 of the slab thick- 
ness. Then the transmission probability through the slab 
depends on i ^ so as 

sq^"-^ a - 1 
T. 

1 - a 



2 — a 
c 



^"'^ 21n(cL/so)' 



2-a 
for a < 1, 

for a — 1 



for 1 < a < 2, 



T — 



gp ln{cL/so) 
2L 



for a 



(Al9a) 
(A19b) 
(Al9c) 

(A19d) 



Hence Teq oc 1/L 



'^(0,01-1) logarithmic corrections 

for a = 1 and a = 2). This is the same scaling as for the 
Levy walk without truncation (see Sec. I A 3p . 



b. NonequiUbrium initial conditions 

For nonequilibrium initial conditions the transition to 
the regular diffusive regime happens while the walker is 
inside the slab. We may therefore assume that the usual 
diffusive scaling Tnoneq — cr/L applies (with a playing 
the role of the mean free path). In view of Eq. (|A13p . an 
anomalous scaling appears if Smax — cL scales propor- 
tionally to L, 



nonoq — (CJ^ ) Sq (X Ij 



(A20) 

The anomalous L^"/^ scaling of Sec. lA 21 now appears as 
a consequence of regular diffusion with a scale dependent 
mean free path. 



T 1 — I I I I I I 




L/so 

FIG. 12: Transmission probability Tnoneq of a Levy walk 
through a slab of thickness L, for nonequilibrium initial con- 
ditions. The data points are the results of a numerical sim- 
ulation, for different values of the step size exponent a (and 
fixed Smax S> L). The lines indicate the expected L~'^^^ scal- 
ing. For a < 1 we also have an analytical prediction (|A5|) for 
the prefactor (solid lines), while for q > 1 only the exponent 
is known analytically so the prefactor has been fitted to the 
data (dotted lines). 
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FIG. 13: Transmission probability Tcq of a Levy walk through 
a slab of thickness L, for equilibrium initial conditions. The 
two panels are for diflPerent values of a. The data points result 
from a numerical simulation, with different values Smax of the 
maximum step size. The solid line is the asymptote ()A9|) for 



Numerical test 



We have tested the analytical expressions (jASP and 
9| by numerical simulation. Results for Tnoneq are 
shown in Fig. 1121 This is the nonequilibrium initial con- 
dition, where the walker starts off at a; = with a step 
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FIG. 14: Transmission probability for a Levy walk with max- 
imum step size Smax that increases proportionally to L. The 
two panels (both for Smax ~ L/10) correspond to equilibrium 
and nonequilibrium initial conditions. The dotted lines show 
the expected scaling (IA19|I and (IA20|I , up to a prefactor which 
has been fitted to the data. (For Tcq the difference with Eq. 
(|A19[I is a factor of two, independent of a.) 



to positive x. The scaling is reproduced for all 

< a < 2, and the prefactor (|A5|) agrees well with the 
simulations for < a < 1. 



For the equilibrium initial condition the walker starts 
off at a large distance from x = 0, crossing the boundary 
at a random point between two scattering events. Re- 
sults of numerical simulations are shown in Fig. 1131 Un- 
like in the nonequilibrium case, the convergence to the 
asymptotic scaling with increasing Smax is very slow, in 
particular for small a. 



We have also tested the scaling (jA19|) and (jA20|) for 
a truncated Levy walk with a maximum step size Smax 
that is a fixed fraction of L. Results are shown in Fig. 
[14] for both equilibrium and nonequilibrium initial condi- 
tions. The anomalous scaling now appears even though 
the diffusion is regular on the scale of L, because of the 
scale dependence of the mean free path. For both types 
of initial conditions the numerics follows closely the ana- 
lytically predicted power laws, including the logarithmic 
factors for a = 1,2 in the equilibrium case. (The con- 
stant prefactors are not given reliably by the analytics.) 
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